Adiabatic cooling of Fermions in an optical lattice 
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The entropy-temperature curves are calculated for non-interacting fermions in a 3D optical lattice. These 
curves facilitate understanding of how adiabatic changes in the lattice depth affect the temperature, and we 
demonstrate regimes where the atomic sample can be significantly heated or cooled. When the Fermi energy of 
the system is near the location of a band gap the cooling regimes disappear for all temperatures and the system 
can only be heated by lattice loading. For samples with greater than one fermion per site we find that lattice 
loading can lead to a large increase in the degeneracy of the system. We study the scaling properties of the 
system in the degenerate regimes and assess the effects of non-adiabatic loading. 
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I. INTRODUCTION 

Tremendous progress has been made in the preparation, 
control and manipulation of Fermi gases in the degenerate 
regime CI S B B S 01 • Such systems have many po- 
tential applications in the controlled study of fermionic su- 
perfluidity and the production of ultra-cold molecules. An- 
other area of developing theoretical interest is in the physics 
of fermions in optical lattices 0, [ijj Hill , and initial ex- 
periments have already begun to examine the properties of 
Fermi-gases (prepared as boson-fermion mixtures) in one- 
dimensional optical lattices II 2l Hal . For Bose gases, opti- 
cal lattices have been used to demonstrate an impressive ar- 
ray of experiments such as: quantum matter-wave engineering 
fl4 [l5ll : the Mott-insulator quantum-phase transition |16]; 
quantum entanglement 1 17]; and coherent molecule produc- 
tion fljil . It seems likely that a similar range of rich physics 
lies ahead for fermions in optical lattices. 

Many of the physical phenomenon that are suitable to ex- 
perimental investigation in optical lattices are sensitive to tem- 
perature and it is therefore of great interest to understand how 
the temperature of a quantum degenerate gas changes with 
lattice depth. Experimental results by Kastberg et al. 11911 
in 1995 showed that loading laser cooled atoms into a three- 
dimensional optical lattice caused the atoms to increase their 
temperature 1 26] . Recently one of us conducted a detailed 
thermodynamic study of bosonic atoms in optical lattices 12011 . 
In that work we showed that for sufficiently low initial tem- 
peratures a new regime would be entered in which adiabati- 
cally ramping up the lattice depth would have the desirable ef- 
fect of cooling the system. The typical temperatures at which 
Bose-Einstein condensates are produced lie well within this 
cooling regime, and thus benefit from reduced thermal fluc- 
tuations when adiabatically loaded into an optical lattice. In 
this paper we examine how degenerate fermions are affected 
by adiabatic loading into an optical lattice. In Fermi gases the 
lowest temperatures obtained in experiments tend to be much 
higher than in Bose gas experiments. It is therefore important 
to understand to what extent the introduction of an optical lat- 
tice might affect the temperature, in particular to determine in 
what regimes additional cooling can occur during lattice load- 
ing. 

The quintessential difference in the properties of degener- 



ate fermions from bosons is embodied by the Fermi energy 
- the energy that marks the top of the Fermi-sea of occupied 
states (at T = 0). The Fermi energy sets a new energy scale 
that has no analog in boson systems and plays a crucial role 
in determining the effect that lattice loading has upon the sys- 
tem. We find that as the Fermi energy approaches a band gap, 
the cooling regime vanishes and the system can only heat with 
increasing lattice depth. However, we also find that when the 
Fermi energy lies in the second band (when the average num- 
ber of fermions per site is greater than one) a cooling regime 
is re-established. This cooling regime for the second band is 
accompanied by a large amplification of degeneracy, i.e. adi- 
abatically loading into the lattice causes both T and the ratio 
T/Tp to decrease. 

The results we present in this paper are obtained from a 
numerical study of the thermodynamic properties of an ideal 
gas of fermions in a 3D cubic lattice. We work with the 
grand canonical ensemble and use the exact single particle 
eigenstates of the lattice to determine the entropy-temperature 
curves for the system for various lattice depths and filling fac- 
tors. We develop analytic expressions for the plateaus that 
develop in the entropy-temperature curves and characterize 
a scaling relationship that holds for low temperatures and in 
deep lattices. A fast loading procedure is considered to ascer- 
tain how robust our results are to non-adiabatic effects. The 
physics we explore here will be relevant to current experi- 
ments, and many of the predictions we make should be easily 



H. FORMALISM 

A. Single Particle Eigenstates 

We consider a cubic 3D optical lattice made from 3 in- 
dependent (i.e. non-interfering) sets of counter-propagating 
laser fields of wavelength A, giving rise to a potential of the 
form 



VLatt(r) 



V 



cos(2fcr) + cos(2ky) + cos(2kz)], (1) 



where k — 2tt/X is the single photon wavevector, and V is the 
lattice depth. We take the lattice to be of finite extent with a 
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Figure 1 : The dependence of the energy gap (e gap dashed line) and 
ground band width (ebw solid line) on the lattice depth (see the text). 



total of N s sites, consisting of an equal number of sites along 
each of the spatial directions with periodic boundary condi- 
tions. The single particle energies e q are determined by solv- 
ing the Schrodinger equation 



e q ^ q (r) = f^q(r) + U Latt (r)</> q (r), (2) 

for the Bloch states, ^> q (r), of the lattice. For notational sim- 
plicity we choose to work in the extended zone scheme where 
q specifies both the quasimomentum and band index of the 
state under consideration [27]. By using the single photon 
recoil energy, Er — h 2 k 2 /2m, as our unit of energy, the en- 
ergy states of the system are completely specified by the lattice 
depth V and the number of lattice sites N s (i.e. in recoil units 
e q is independent of k). 

For completeness we briefly review some important fea- 
tures of the band structure of Eq. (0 relevant to the thermody- 
namic properties of the system. For sufficiently deep lattices 
an energy gap, e gap , will separate the ground and first excited 
bands (see Fig. [Q. For the cubic lattice we consider here, a 
finite gap appears at a lattice depth of V « 2Er t2M (marked 
by the vertical asymptote of the dashed line in Fig. [Q. For 
lattice depths greater than this, the gap increases with lattice 
depth. In forming the gap higher energy bands are shifted up- 
wards in energy, and the ground band becomes compressed — 
a feature characteristic of the reduced tunneling between lat- 
tice sites. We refer to the energy range over which the ground 
band extends as the (ground) band width, 6bw- As is apparent 
in Fig. ^ the ground band width decreases exponentially with 
V , causing the ground band to have an extremely high density 
of states for deep lattices. 



B. Equilibrium Properties 

Our primary interest lies in understanding the process of 
adiabatically loading a system of N p fermions into a lattice. 



Under the assumption of adiabaticity the entropy remains con- 
stant throughout this process and the most useful information 
can be obtained from knowing how the entropy depends on 
the other parameters of the system. In the thermodynamic 
limit, where N s — > oo and N p — ■> oo while the filling fac- 
tor n = N p /N s remains constant, the entropy per particle 
is completely specified by the intensive parameters T, V, n. 
The calculations we present in this paper are for finite size 
systems, that are sufficiently large to approximate the thermo- 
dynamic limit. We would like to emphasize at this point the 
remarkable fact that V is an adjustable parameter in optical 
lattice experiments, in contrast to solid state systems where 
the lattice parameters are immutable. 

The entropy is determined as follows: The single particle 
spectrum {e q } of the lattice is calculated for given values of 
N s and V. We then determine the thermodynamic properties 
of the lattice with N p fermions in the grand canonical ensem- 
ble, for which we calculate the partition function Z 

logZ = ^log(l + e-^-")), (3) 
q 

where /i is found by ensuring particle conservation. The en- 
tropy of the system can then be expressed as 

S = k B (log Z + (3E - npN p ) , (4) 

where f3 = 1 /fcsT, and E = —d In Z/d/3 is the mean energy. 

1. Multiple components 

In most current experiments mixtures of Fermi gases in dif- 
ferent internal states are studied. This is required because s- 
wave elastic collisions, needed for re-equilibration, are pro- 
hibited by the Pauli principle for spin polarized samples 1 29] . 
The theory we present here is for the spin polarized case, but is 
trivially extensible to multiple components if the lattice poten- 
tial is spin independent and the number of atoms in each com- 
ponent is the same: in this case all extensive parameters are 
doubled (e.g. {E, S}) and intensive parameters (e.g. {T, /i}) 
remain the same. The inclusion of interaction effects, which 
will be important in the multiple component case, is beyond 
the scope of this paper. 



III. RESULTS 
A. Effect of lattice loading on Fermi-gas temperature 

In Fig. |2]we show entropy-temperature curves for various 
lattice depths and filling factors n. These curves have been 
calculated for a lattice with 31 lattice sites along each spatial 
dimension, i.e. N s w 3 X 10 4 . 

A general feature of these curves is the distinct separation 
of regions where adiabatic loading causes the temperature of 
the sample to increase or decrease, which we will refer to as 
the regions of heating and cooling respectively. These regions 
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Figure 2: Entropy versus temperature curves for a N s « 3 x 10 4 site cubic lattice, at various depths V — to 20Er (with a spacing of 2Er 
between each curve). Filling factors used are (a) n — 0.25, (b) n — 0.8, (c) n — 1.0, and (d) n = 1.2. The entropy plateau is shown as a 
dashed line. Dotted line marked A shows a path along which adiabatic loading into the lattice causes the temperature to increase. Dotted line 
marked B shows a path along which adiabatic loading into the lattice causes the temperature to decrease. 



are separated by a value of entropy at which the curves plateau 
- a feature that is more prominent in the curves for larger lat- 
tice depths. This plateau entropy is indicated by a horizontal 
dashed line and is discussed below. For the case of unit filling 
factor shown in Fig. EJc), this plateau occurs at S = 0, and 
only a heating region is observed. 

We now explicitly demonstrate the temperature changes 
that occur during adiabatic loading using two possible adia- 
batic processes labeled A and B, and marked as dashed-dot 
lines in Fig. |2ja). Process A begins with a gas of free parti- 
cles in a state with an entropy value lying above the plateau 
entropy. As the gas is loaded into the lattice the process line 
indicates that the temperature increases rapidly with the lattice 
depth. Conversely process B begins with a gas of free parti- 
cles in a state with entropy below the plateau. For this case 



adiabatic lattice loading causes a rapid decrease in tempera- 
ture. This behavior can be qualitatively understood in terms of 
the modifications the lattice makes to the energy states of the 
system. As is apparent in Fig. ^ the ground band rapidly flat- 
tens for increasing lattice depth causing the density of states 
to be more densely compressed at lower energies. Thus in the 
lattice all these states can be occupied at a much lower temper- 
ature than for the free particle case. As we discuss below, So 
is the maximum entropy available from only accessing states 
of the lowest band. If S < So, the temperature of the sys- 
tem must decrease with increasing lattice depth to remain at 
constant entropy. Alternatively, for S > So the occupation of 
states in higher bands is important, and as the lattice depth and 
hence e gap increases, the temperature must increase for these 
excited states to remain accessible. 
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B. Fermi-gas degeneracy 

In addition to the effect that lattice loading has on the tem- 
perature of a Fermi-gas, it is of considerable interest to under- 
stand how the ratio of temperature to the Fermi temperature 
(Tf)|30[ changes. Indeed, the ratio T/Tp is the standard fig- 
ure of merit used to quantify the degeneracy of dilute Fermi 
gases. In Fig. [3]we show how T/Tf changes with adiabatic 
lattice loading for the same parameters used in Fig. [3 In Fig. 
I3(a)-(b) the same general behavior is seen: Below the entropy 
plateau where cooling is observed (see Fig. |2(a)-(b)), the ra- 
tio of T/Tp remains approximately constant, so that there is 
little change in the degeneracy of the gas. Above the entropy 
plateau where heating was observed, the ratio of T/Tp rapidly 
increases, so that in this regime the gas will rapidly become 
non-degenerate as it is loaded into the lattice. For the unit 
filling case (Fig. 13(c)), there is no cooling regime, and heat- 
ing is accompanied by a rapid increase in T/Tp for all initial 
conditions of the gas. In Fig. 0(d), where the filling factor 
is n = 1.2, rather different behavior is seen: In the cooling 
regime, the ratio of T/Tp is rapidly suppressed as the temper- 
ature decreases, e.g. see the dotted line marked C in Fig. [5] 
(d). This most desirable behavior could be used for example 
to prepare a Fermi-gas into a highly degenerate state where 
the BCS transition might be observable. We also note that for 
the same parameters but in the heating regime the ratio T/Tp 
remains relatively constant. 

We can give a simple explanation for the behavior of T/Tp. 
For the three cases considered in Fig. |3](a)-(c) the Fermi en- 
ergy lies within or at the top of the first band of energy states. 
As shown in Fig. ^ me width of the ground band (eew) de- 
creases rapidly with lattice depth. Because the number of 
states contained in each band is constant (given by the number 
of lattice sites) both the Fermi energy and Tp scale identically 
to ebw, thus will rapidly decrease with lattice depth. In the 
cooling regime, the temperature scales in the same manner as 
£bw ( see Sec. mi Dl below and Fig. |4ji, and thus the ratio 
T/Tp remains approximately constant. In the heating regime 
T increases slowly, while the ratio T/Tp increases rapidly 
with lattice depth (on account of Tp becoming small). 

For the case considered in Fig. |5](d) the filling factor satis- 
fies n > 1 and the Fermi energy lies in the second band. As 
the lattice depth increases the Fermi energy and Tp now scale 
like e gap , i.e. slowly increases with lattice depth (see Fig. 0. 
Thus in the regime where the temperature decreases, the ratio 
T/Tp must become smaller. We note that the temperature re- 
duction occurs because the width of the second band decreases 
with lattice depth. 



C. Entropy plateau 

In Figs. |2ja) and (b) a horizontal plateau (at the level 
marked by the dashed lines) is common to the entropy- 
temperature curves for larger lattice depths [V > 8Er). This 
occurs because for these lattices, the energy range over which 
the ground band extends is small compared to the energy gap 
to the excited band, and there is a large temperature range 



over which states in the excited bands are unaccessible, yet 
all the ground band states are uniformly occupied. The en- 
tropy value indicated by the dashed line in Figs. [2a) and (b) 
corresponds to the total number of ATp -particle states in the 
ground band. Since the number of single particle energy states 
in the ground band is equal to the number of lattice sites, the 
total number of available A^-particle states (Slo) is given by 
fio = N s \/[N p \(N s - N p )\] (valid for N p < N s ). The asso- 
ciated entropy So = fcs log £7q, which we shall refer to as the 
plateau entropy, can be evaluated using Sterling's approxima- 
tion 

So =s k B [N p logN p + N s logN s (5) 
-(2V s -Agiog (N s -N p ) 

the validity condition for this result is that 1 <C N s -C N p . An 
important case for which the above approximation is invalid 
is for N p — N B , i.e. we have a filling factor of n = 1, where 
So = 0. This case corresponds to the unit filling factor result 
shown in Fig. |3c) where, as a result of the entropy plateau 
occurring at S = 0, only a heating region is observed. 

Similar entropy plateaus are observed for greater than unit 
filling (N p > N s ), e.g. as is seen in Fig. |2jd). For fermions 
such high filling factors necessarily means that higher bands 
are occupied, and in general the precise details of these higher 
plateaus will depend on the particle band structure of the lat- 
tice. E.g in the lattice we consider here Q there are three 
degenerate first excited bands, that contain a total of 3N S sin- 
gle particle states. Because the first band is fully occupied 
only N p — N p — N s particles are available to occupy the ex- 
cited band, and so the total number of available states is found 
according to the ground band result (|5jl, but with the substitu- 
tions N s — > 3N S and N p — > N p — N s . This result is shown as 
the dashed horizontal in Fig. [3d) labelled as Si. 

The suppression of the plateaus at specific integer filling 
factors (e.g. n = 1 for So — > and n = 4 for Si — > 0) corre- 
sponds to the Fermi-energy of the system approaching a band 
gap. Whenever this occurs it means that all the states below 
the gap are occupied at T = and excitations in the system 
require the promotion of particles into the excited band (above 
the gap). As all band gaps increase in size with lattice depth, 
the temperature of the system must increase for the entropy to 
remain constant. Thus in regimes where the Fermi-energy lies 
at a band gap the system only exhibits heating with increasing 
lattice depth (e.g. see Fig. 13c)). 



D. Scaling: Tight-binding limit at low temperatures and filling 
factors 



Here we give limiting results for the entropy-temperature 



As discussed in Sec. MI CI when N p < N s and the temper- 
ature is sufficiently low that S < Sq, then only single particle 
states within the ground band are accessible to the system. 
In addition when the tight-binding description is applicable 
for the initial and final states of an adiabatic process, the ini- 
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Figure 3: Entropy versus temperature over Fermi temperature curves for the same cases considered in Fig. |5|(see that figure for details of 
parameters used). The entropy plateau is shown as a dashed line. The dotted line marked C shows a path along which adiabatic loading into 
the lattice causes the ratio of the temperature to the Fermi temperature to decrease. 



tial and final thermodynamic variables are related by a scaling 
transformation. 

In the tight-binding regime, which is good approximation 
f° r V ^Er, the ground band dispersion relation takes the 
form 

£TB(q) = -^ cos (* a )> ( 6 ) 

j={x,y,z} 

where a = 7r/A is the lattice period, the ground bandwidth 
£bw has already been introduced (e.g. see Fig. 0, and the 
wavevector q is restricted to the first Brillouin zone. We re- 
fer the reader to Refs. l2ll E2I1 for more details on the tight- 
binding approximation. 

To illustrate the scaling transformation we consider an ini- 
tial system in equilibrium with entropy S < So, in lattice of 
depth Vi sufficiently large enough for tight-binding expression 



to provide an accurate description of the ground band en- 
ergy states. If an adiabatic process is used to take the system 
to some final state at lattice depth Vf (also in the tight-binding 
regime) it is easily shown that the macroscopic parameters of 
the initial and final states are related as 

X f = aX u (7) 

where X = {E, T, or /i}, and the scaling parameter a = 
(eBw)//(cB\v)i is given by the ratio of the final and initial 
bandwidths. The requirement that the initial and final states 
are in the tight-binding regime is because the single particle 
states are then related as (era (l) ) / = ct (era(q))j, which is 
essential for to hold. 

This type of scaling suggests that the occupations of the sin- 
gle particle levels are unchanged during the change in lattice: 
the products /3exB(q) and (3[i are independent of V, so the 



6 



10' 



10 



DQ o 

£ 10 



10" 



10" 



T.=0.45, S/S =1.1 

i o 



T.=0.35, S/S =0.97 
i o 



T.=0.2, S/S =0.66 
i o 



T.=0.05, S/S =0.19 
i o 



x 10 



5 10 15 

V [E R ] 



20 



Figure 4: Scaled temperature change during adiabatic lattice load- 
ing. Ratio of temperature to bandwidth for various lattice depths 
along adiabatic contours for several values of entropy. Initial tem- 
perature (at V = 0) and ratio of entropy to the plateau entropy is 
indicated for each curve. Results are indicated for the case of filling 
factor n = 0.25. 



Fermi distribution, /f(q) = [exp(/3exB(q) — /3/J,) + will 
also be independent of V. This suggests that being adiabatic 
in this regime will not require redistribution through collisions 
and may allow the lattice depth to be changed more rapidly. 

To confirm the scaling predictions, in Fig. @]we plot the 
ratio of the temperature to ground bandwidth as a function of 
lattice depth along contours of constant entropy, (e.g. how 
/cbT/ebw varies along the process curves labelled A and B in 
Fig. EJa)). In regimes where the scaling relationship Q holds 
true, the ratio /cbT/cbw should be constant (independent of 
V). In Fig. |4]this is clearly observed for initial entropies less 
than So and lattice depths V > 5-Er. For S > Sq single 
particle states of higher bands necessarily play an important 
role in the thermodynamic state of the system, and the scaling 
transformation clearly does not hold at any lattice depth, as 
is seen in the dotted curve in Fig. |4] For this case as the 
lattice depth increases the cooling effect of the ground band 
compression is offset by the f particles in the excited band 
that are lifted to larger energies as the gap (e gap ) grows (see 
Fig-E). 




12 3 

Temperature: ^T/E R 

x 10 4 




Temperature: i^T/E 



'R 



Figure 5: Fast lattice loading of a N s ~ 3 x 10 4 site cubic lat- 
tice, with filling factors of (a) n — 0.25 and (b) n = 1. Dark solid 
lines indicate fast loading curves (see text). The initial V — state 
for each of this curves is indicated with a dot. The lattice depth on 
these curves can be determined from their intercept with the equilib- 
rium entropy versus temperature curves (gray solid lines), which are 
described in Figs. 0a) and (c) for n = 0.25 and n = 1 respectively. 



E. Adiabaticity 

Finally we note that interactions between particles are es- 
sential for establishing equilibrium in the system, and un- 
derstanding this in detail will be necessary to determine the 
timescale for adiabatic loading. In general this requirement 
is difficult to assess, and in systems where there is an addi- 
tional external potential it seems that the adiabaticity require- 
ments will likely be dominated by the process of atom trans- 
port within the lattice to keep the chemical potential uniform, 
though recent proposals have suggested ways of reducing this 



problem l23ll for Bose systems. A study of the effects of in- 
teractions or inhomogeneous potentials is beyond the scope of 
this work, however it is useful to assess the degree to which 
non-adiabatic loading would cause heating in the system. We 
consider lattice loading on a time scale fast compared to the 
typical collision time between atoms, yet slow enough to be 
quantum mechanically adiabatic with respect to the single par- 
ticle states. This latter requirement excludes changing the lat- 
tice so fast that band excitations are induced, and it has been 
shown that in practice this condition can be satisfied on very 
short time scales 1 24] . We will refer to this type of loading as 
fast lattice loading, to distinguish it from the fully adiabatic 
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loading we have been considering thus far. 

To simulate the fast lattice loading we take the system 
to be initially in equilibrium at temperature T{ for zero lat- 
tice depth. For the final lattice depth we fast load into, we 
map the initial single particle distribution onto their equiva- 
lent states in the final lattice, and calculate the total energy 
for this final non-equilibrium configuration (i.e. we calculate 

E = E q £ q (/) Meq\ Zi), where e q « and e q ^) are the sin- 
gle particle energies for the initial and final lattice depths re- 
spectively, and fp is the Fermi distribution function). This 
procedure assumes that there has been no collisional redistri- 
bution to allow the system to adjust to the lattice potential dur- 
ing the period it is changed. To determine the thermodynamic 
state the final distribution will relax to, we use the energy of 
the non-equilibrium distribution as a constraint for finding the 
equilibrium values of temperature and entropy. In general the 
final state properties will depend on the initial temperature, 
filling factor, and final depth of the lattice. To illustrate typ- 
ical behavior we show a set of fast loading process curves in 
Fig. |5]for two different values of the filling factor. 

These curves show, as is expected from standard thermo- 
dynamic arguments, that entropy increases for non-adiabatic 
processes, i.e. all loading curves in Figs. 0a) and (b) bend up- 
wards with increasing lattice depth. For the results with filling 
factor n = 0.25 and for initial temperatures deep in the cool- 
ing regime (i.e. initial states far below the entropy plateau) a 
useful degree of temperature reduction can be achieved with 
fast lattice loading up to certain maximum depth. For exam- 
ple, the second lowest fast loading curve in Fig. |3Ja) cools 
with increasing lattice depth up to V ~ 15Er, and then be- 
gins to heat for larger final lattice depths. Generally for low 
filling factors (n < 1) where the ground band plays the dom- 
inant role in the system behavior at low temperatures, the en- 
tropy increase is due to mainly to the reshaping of the single 
particle energy states which occurs at low lattice depth [31]. 
This effect can be reduced by taking as the initial condition 
for the fast loading a system in equilibrium at a finite lattice 
depth for which the dispersion relation is more tight-binding- 
like. This situation was considered in Ref. |8[ in application 
to preparing a superfluid Fermi gas using in an optical lattice. 
Their results, for the case n — 0.5 and an initial lattice depth 



of V w IEr, predicted a useful degree of cooling. 

As was demonstrated in Fig. |2jc), for filling factor n = 1 
adiabatic lattice loading causes the atoms to heat. This effect 
is exasperated by non-adiabatic loading, as shown in Fig. EJb). 
This case also benefits from beginning in a lattice of non-zero 
depth, since at fixed temperature but increasing lattice depth 
(hence larger €bw), a smaller number of particles will be found 
in the excited bands. 



IV. CONCLUSION 

In this paper we have calculated the entropy-temperature 
curves for fermions in a 3D optical lattice at various depths 
and filling factors. We have identified general features of 
the thermodynamic properties relevant to lattice loading, indi- 
cated regimes where adiabatically changing the lattice depth 
will cause heating or cooling of the atomic sample, and 
have provided limiting results for the behavior of the en- 
tropy curves. The results presented in this work suggest op- 
timal regimes (filling factors and temperatures) which will fa- 
cilitate the suppression of thermal fluctuations in a fermionic 
gas by lattice loading. These predictions should be easily ver- 
ifiable with current experiments. We have also shown that 
for a sample of fermions with filling factor greater than one, 
the cooling regime is accompanied by a significant reduction 
of the temperature compared to the Fermi temperature. This 
regime would clearly be desirable for experiments to inves- 
tigate as an avenue for producing dilute Fermi gases with 
T/Tp <C 1. We have shown that many of our predictions 
are robust to non-adiabatic effects. 
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